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Abstract 



In this paper we dynamically determine the mass of the Kuiper 
Belt Objects by exploiting the latest observational determinations of 
the orbital motions of the inner planets of the Solar System. Our result, 
in units of terrestrial masses, is 0.033±0.115 by modelling the Classical 
Kuiper Belt Objects as an ecliptic ring of finite thickness. A two-rings 
model yields for the Resonant Kuiper Belt Objects a value of 0.018 ± 
0.063. Such figures are consistent with recent determinations obtained 
with ground and space-based optical techniques. Some implications 
for precise tests of Einsteinian and post-Einsteinian gravity are briefly 
discussed. 
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1 Introduction 

Starting in 1992, astronomers have become aware of a vast population of 
small bodies orbiting the Sun beyond Neptun^|. There are at least 70,000 
Trans-Neptunian Objects (TNOs) with diameters larger than 100 km in the 
30-50 AU region. Observations show that TNOs are mostly confined within 
a thick band around the ecliptic, leading to the realization that they occupy 
a ring or belt surrounding the Sun. This ring is generally referred to as the 
Kuipei@Belt (Edgeworth 1943; Kuiper 1951; Fernandez 1980). 
Reasons of interest in the Kuiper Belt are as follows 



1 See http://www.ifa.hawaii.edu/faculty/jewitt/kb.html] 

2 See http://www.ifa.hawaii.edu/faculty/jewitt/kb/gerard.html for the origin of the 
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• It is likely that the Kuiper Belt Objects (KBOs) are extremely prim- 
itive remnants from the early accretional phases of the solar system. 
The inner, dense parts of the pre-planetary disk condensed into the ma- 
jor planets, probably within a few millions to tens of millions of years. 
The outer parts were less dense, and accretion progressed slowly. Ev- 
idently, many small objects were formed 

• It is widely believed that the Kuiper Belt is the source of the short- 
period comets. It acts as a reservoir for these bodies in the same way 
that the Oort Cloud acts as a reservoir for the long-period comets 

• KBOs are usually not yet modelled in the data reduction softwares 
used for orbit determination purposes, so that they may represent a 
serious bias in precise tests of Einsteinian and post-Einsteinian gravity. 
Thus, it is important to assess their impact on planetary motions 

KBOs can be classified into three dynamical classes (Jewitt et al. 1998) 

• Classical KBOs (CKBOs), following nearly circular orbits with rela- 
tively low eccentricities (e < 0.25) and semimajor axes 41 AU < a < 
46 AU; they constitute about 70% of the observed population 

• Resonant KBOs occupy mean motion resonances with Neptune, such 
as 3:2 (the Plutinos, a ~ 39.4 AU) and 2:1 (a ~ 47.8 AU), and amount 
to about 20% of the known objects 

• Scattered KBOs, which represent only about 10% of the known KBOs 
but have the most extreme orbits, with a ~ 90 AU and e ~ 0.6, 
presumably due to a weak interaction with Neptune; we only have 
rather poor knowledge of them. 

In this paper we determine the mass of KBOs in a truly dynamical way by 
means of the latest observational determinations of the secular perihelion 
advances of the inner planets of the Solar System. The use of such celes- 
tial bodies presents many advantages. Long data sets including also many 
accurate radio-technical range and range-rate measurements are available; 
they cover many orbital revolutions, contrary to the outer planets which are 
more affected by KBOs like Uranus and Neptune. Indeed, mainly optical 
data and sparse radar-ranging measurements exist for them covering barely 
one period for Uranus and less than one full orbital revolution for Neptune. 
Moreover, for heliocentric distances of the order of just 1 AU or less many 
details of the true mass distribution of KBOs are not relevant and all the re- 
cently proposed analytical models converge satisfactorily to a substantially 
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Table 1: Observationally determined extra-precessions of the longitudes of 
perihelia of the inner planets, in arcseconds per century (" cy" 1 ), by using 
EPM2004 with /3 = 7 = 1, J 2 = 2x 1(T 7 , from Table 3 of (Pitjeva 2005b). 
Both KBOs and the general relativistic gravitomagnetic force were not in- 
cluded in the adopted dynamical force models. The quoted uncertainties are 
not the mere formal, statistical errors but are realistic in the sense that they 
were obtained from comparison of many different solutions with different 
sets of parameters and observations. The correlations among such deter- 
mined planetary perihelia rates are very low with a maximum of about 20% 
between Mercury and the Earth (Pitjeva, private communication 2005). 
Mercury Venus Earth Mars 

-0.0036 ± 0.0050 0.53 ±0.30 -0.0002 ± 0.0004 0.0001 ± 0.0005 



unified description, given the present-day accuracy in reconstructing the or- 
bital motions in our region of the Solar System. This fact greatly simplifies 
the analytical calculation and allows for reliable, rather model-independent 
determinations of the KBOs's mass. 

2 The determined extra-rates of perihelion 

The Russian astronomer E.V. Pitjeva (Institute of Applied Astronomy, Rus- 
sian Academy of Sciences) recently processed almost one century of data 
of all types in the effort of continuously improving the EPM2004 plane- 
tary ephemerides (Pitjeva 2005a). Among other things, she also determined 
residual secular, i.e. averaged over one orbital revolution, rates of the longi- 
tudes of perihelion w = Q,+u, where u and ft, are the argument of perihelion 
and the longitude of the ascending node, respectively, of the inner planets 
(Pitjeva 2005b) as fit-for parameters of global solutions in which she con- 
trasted, in a least-square way, the observations (ranges, range-rates, angles 
like right ascension a and declination S, etc.) to their predicted values com- 
puted with a complete suite of dynamical force models. The results are 
shown in Table [H The modelled forces are 

• The Newtonian N-body perturbations including also the effect of 301 
largest asteroids (Krasinsky et al. 2002) and of the minor asteroid ring 
in the ecliptic 

• The Sun's quadrupole mass moment J2 (Paterno et al. 1996; Pijpers 
1998; Mecheri et al. 2004), set to 2 x 10~ 7 
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• The post-Newtonian gravitoelectric forces (Newhall et al. 1983) pa- 
rameterized in terms of the PPN Eddington-Robertson-Schiff param- 
eters 7 and (3 (Will 1993) set to their general relativistic values 7 = 
(3 = 1 

The un-modelled forces are 

• The Newtonian force induced by KBOs 

• The post-Newtonian gravitomagnetic force responsible for the Lense- 
Thirring effect (Lense and Thirring 1918). 

Thus, the effect of KBOs is entirely accounted for by the so-obtained residual 
perihelia advances of Table [U 



3 Modelling and confrontation with data 

In order to make a comparison with Table [TJ a theoretical prediction for the 
secular precessions induced by KBOs on the planetary perihelia is required. 

The action of KBOs can be treated in a perturbative way. The Gauss 
rate equations for lo and Q are, for a generic orbital configuration 



~dt 



dtt Vl - ( 

cos i— — I 

at nae 
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(1) 



and 



~dt 
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where i is the the inclination of the orbit, p = a(l — e 2 ), n = \J GM/a 3 is the 
Keplerian mean motion, / is the true anomaly counted from the pericentre, 
and A r , At and A n are the radial, transverse and normal components of the 
perturbing acceleration A, respectively. 

For almost ecliptic orbits {i ~ deg, cosi ~ 1, sini ~ i), like those 
of the Solar System's major bodies, the rate equation for w can safely be 
approximated as 



dm 
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Its secular rate is obtained by evaluating the right-hand-side of eq. (J3J) onto 
the unperturbed Keplerian ellipse, characterized in terms of the eccentric 
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anomaly E by 



r = a(l — e cos E), 

«*»/ = £Sss. (4) 



and averaging the result over one orbital period P by means of 

dt = (1-ecosE) 

P 2vr W 

Motivated by the search for a gravitational explanation of the Pioneer 
anomaly (Anderson et al. 2002), various authors recently produced several 
analytical models for the KBOs acceleration (Anderson et al. 2002; Nieto 
2005; Bertolami and Vieira 2006; de Diego et al. 2006). At heliocentric 
distances < 1 AU many details of the three-dimensional mass distribution of 
KBOs can be neglected, so that we can approximate it with a bi-dimensional 
distribution lying in the ecliptic plane. Moreover, the four models analyzed 
by Bertolami and Vieira (2006), i.e. two-rings, uniform disk, non-uniform 
disk and torus, give the same results for r <C 20 AU and for various ecliptic 
latitudes (3 close to zero. Thus, we will adopt for CKBOs the uniform disk 
model, which consists of a uniform, hollow thin disk lying in the ecliptic 
plane within distances i2 m i n = 37.8 AU and -R max = 46.2 AU (Bernstein et 
al. 2004). 

The Newtonian gravitational field of a thin massive ring exhibits cylin- 
drical symmetry and has an in-plane, radial component and an out-of-plane, 
normal component (Owen 2003). Since we are looking for the effects induced 
on the longitude of the perihelion, only the radial component is relevant to 
us. Thus, from eq. ([3])— eq. © the secular perihelion rate can be written as 



dt / 2vre V GM 



dw\ 1 /a(l-e 2 ) f 2n 



/ A r (E)(cosE-e)dE. (6) 
Jo 



The radial component of Ak for a uniform disk is (Bertolami and Vieira 
2006) 



2tt r R 
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Table 2: Gravitomagnetic secular precessions of the longitudes of perihelion 
w of Mercury, Venus, Earth and Mars in " cy" 1 . The value (190.0 ± 1.5) x 
10 39 kg m 2 s" 1 (Pijpers 1998; 2003) has been adopted for the solar proper 
angular momentum Lq. 

Mercury Venus Earth Mars 

-0.0020 -0.0003 -0.0001 -0.00003 



rriK is the KBOs's mass, i?x an d 4>k are the polar coordinates of the disk 
mass element, and {r, /3, A} are the usual spherical ecliptic planetary coor- 
dinates. For the inner planets, with r <C Rmim -Rmax and whose orbits lie 
almost exactly in the ecliptic plane, we can safely pos^l 
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yk) - 2 (^)cos/?cos(0k-A) 



- 3 / 2 ( v 

~ 1 + 3 i^-j cos(0k-A). 



(9) 

By inserting eq. Q into eq. (|7|) and performing the integration we get 

Gm-u 

A nd = — r. (10) 

2(-Rmax ~i~ ^minj^max-^min 

Note that it is positive, in according with the behavior of all the models 
considered in (Bertolami and Vieira 2006) for r <C 20 AU. Eq. (110j) . with 
r = a(l — ecos-E), into eq. © finally yields 



dw\ _ 3 /Ga 3 (l - e 2 ) m K 



dt I 4 V M (-Rmax -Rmin)-Rmax-^mi]i 



Table [T] and eq. (jlip can now be used to determine dik in a truly 
dynamically, model-independent way. To this aim, let us note that, by 
construction, the determined extra-rates of Table Q] are not only due to 
KBOs, but also the Lense-Thirring field and the mismodelled part of the 
solar J2, which is presently uncertain at a ~ 10% level, contribute to them. 
Their nominal magnitudes are listed in Table [2] and Table respectively. 

Although the Lense-Thirring effect and the mismodelled part of the J2 
precessions may be neglected in our analysis, being smaller than the errors of 
Table [U we prefer conservatively to cancel out, by construction, any possible 
impact due to them. It can be done by suitably combining the perihelia of 
Mercury, the Earth and Mars according to an approach followed, e.g., in 



3 Recall that cos f3 = yT— sin 2 i sin 2 /. 
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Table 3: Nominal values of the classical secular precessions of the longitudes 
of perihelion w of Mercury, Venus, Earth and Mars, in " cy" 1 , induced by 
the solar quadrupolar mass moment J 2- The value J2 = 2 x 1CP 7 used in 
(Pitjeva 2005b) has been adopted. Their mismodelled amplitudes can be 
obtained by assuming an uncertainty in J2 of the order of ~ 10%. 

Mercury Venus Earth Mars 

0.0254 0.0026 0.0008 0.0002 



(Iorio 2005). By converting Tabled] in s , we have the CKBOs' mass, in 
units of terrestrial masses 

(comb) _ ^Mercury + CjCTEarth + C2?Z7Mars _ n nr -r, s 

™ K " 1.017 x 10-15 s -i - U - U52 ' 

with 

d = -81.71, c 2 = 221.18. (13) 

Such coefficients, built in terms of the semimajor axes and eccentricities of 
the involved planets, assure that the solar Newtonian quadrupolar and post- 
Newtonian gravitomagnetic fields, whatever their contributions to Table [1] 
is, do not affect at all the recovered CKBOs' mass. Because of the existing 
correlations among the determined extra-rates of perihelia, the error can be 
conservatively evaluated as 

<■ (comb) , ^Mercury + [cij^gEarth + [cgljgMag _ n 99 o f , 

dm K < 1.017 x 10- 15 s- 1 ~ { ' 

If we only use the Mars perihelion, whose Lense-Thirring and mismod- 
elled J 2 precessions are negligible by one order of magnitude, we get 

m J^ ars) = 0.026 ±0.134. (15) 

A weighted mean of such two measurements yields 

m (weighted) = 033 ± Q n5 

In regard to the Resonant KBOs, we can adopt the two-rings model with 
Ri = 39.4 AU and R2 = 49.8 AU. The radial acceleration is (Bertolami and 
Vieira 2006) 

A 2R - GmK — ^ - Ri COS - cos ^ K ~ ^ dcpK 

2tt(R 1 + R 2 ) J [ r 2 + R 2 _ 2rR . 

cos/3cos(</>k — A)] 

(17) 
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By using the same approximation of eq. ([9]), we obtain mj^° m = 0.029 db 
0.123 and mjf ars) = 0.015 ± 0.073, with m (^ ei s hted ) = 0.018 ± 0.063. Note 
that such figures are smaller than those for CKBOs, consistently with the 
existing estimates for the population of Resonant KBOs. 

At this point, we can a posteriori justify the use of the approximation 
of eq. ([9]) for both the models used: indeed, by using our values for uik it 
is possible to show that additional terms in the expansion of eq. ([9]) yield 
precessions far too small to be detected with the present-day accuracy. 

4 Discussion and conclusions 

In this paper we dynamically determined the mass of KBOs from an analysis 
of the recently determined secular perihelion advances of the rocky planets 
of the Solar System. We modelled CKBOs as a uniform thick ring; such a 
simple model is justified by the fact that at heliocentric distances of about 
1 AU many details of the true three-dimensional KBOs mass distribution 
can be neglected; indeed, many different, more or less complicated analytical 
models manifest substantially the same behavior at distances much less than 
20 AU. For CKBOs we obtained a mass of 0.033±0.115, in units of terrestrial 
masses. A two-rings model for the Resonant KBOs yields a mass of 0.018 ± 
0.063. 

Such figures are consistent with those by 

• Bernstein et al. (2004), who give a nominal mass of 0.010 for CKBOs 
and 0.021 for their Excited class including some Plutinos and Scattered 
KBOs. They used the ACS camera of the Hubble Space Telescope 

• Gladman et al. (2001) yielding a mass of 0.04 - 0.1 for all the TNOs 
in the range 30-50 AU, excluding the Scattered KBOs. The adopted 
technique was deep imaging on the Canada-France-Hawaii Telescope 
and the ESO Very Large Telescope UT1 

• Trujillo et al. (2001) giving for CKBOs a value of 0.030, from a 
wide-field survey with the CCD Mosaic camera of the Canada-France- 
Hawaii Telescope. 

The upper limit of 0.3 terrestrial masses obtained by Backman et al. (1995) 
from far-IR emission measurements is, instead, ruled out. 

It is important to note that all the previous estimates are based on 
various assumptions about, e.g., albedo and density of the KBOs for which 
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Table 4: Nominal values of the classical secular precessions of the longitudes 
of perihelion w of Mercury, Venus, Earth and Mars, in " cy" 1 , induced by 
CKBOs. For wik the value of eq. (fT2|) has been used. 
Mercury Venus Earth Mars 

0.00002 ± 0.0001 0.00006 ± 0.0003 0.0001 ± 0.0004 0.0002 ± 0.0008 



very large uncertainties still exist, so that the authors of the previously cited 
works decided to release no errors of their mass estimations. 

The impact of KBOs on the inner planets of the Solar System lies at 
the edge of the present-day accuracy. Its induced effects, shown in Table H] 
for a particular value of the CKBOs' mass, should not be neglected in, e.g., 
precision tests of gravity because, especially for the Earth and Mars, they are 
of the same order of magnitude, or even larger, than some Einsteinian (see 
Table EJ and post-Einsteinian features of motion^ which recently attracted 
much attention in view of a possible detection in the near future. If not 
accounted for in the dynamical force models of the orbit data reduction 
softwares, KBOs may bias the recovery of such effects when the required 
precision level will be finally attained. In the case of the measurement of the 
Lense-Thirring effect by only using the perihelion of Mercury (Iorio 2005a) , 
the impact of KBOs is negligible. It is not so if a combination of the perihelia 
of Mercury and the Earth is used in order to cancel out the mismodelling 
of J2 (Iorio 2005a). Thanks to the availability of a simple and reliable 
formula as that of eq. (jlip for the inner planets, such a problem could be 
circumvented by setting a suitable three-elements combination allowing for 
de-coupling the Lense-Thirring effect from J2 and KBOs as well. 

Incidentally, let us note that our results further enforce the conclusion 
that the Pioneer anomaly cannot be due to KBOs (Anderson et al. 2002; 
Nieto 2005; Bertolami and Vieira 2006). 
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